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Abstract 

Alms. HESS J 1632-478 is an extended and still unidentified TeV source in the galactic plane. 

Methods. In order to identify the source of the very high energy emission and to constrain its spectral energy distribution, 

we used a deep observation of the field obtained with XMM-Newton together with data from Molonglo, Spitzer and Fermi to 

detect counterparts at other wavelengths. 

Results. The flux density emitted by HESS J 1632-478 peaks at very high energies and is more than 20 times weaker at all other 

wavelengths probed. The source spectrum features two large prominent bumps with the synchrotron emission peaking in the 

ultraviolet and the external inverse Compton emission peaking in the TeV. HESS J 1632-478 is an energetic pulsar wind nebula 

with an age of the order of 10^ years. Its bolometric (mostly GeV-TeV) luminosity reaches 10% of the current pulsar spin down 

power. The synchrotron nebula has a size of 1 pc and contains an unresolved point-like X-ray source, probably the pulsar with 

its wind termination shock. 
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1. Introduction 

HESS J1632-478 has been discovered as a diffuse very high 
energy (VHE) 7-ray source during the 2004-2006 survey 
of the inner Galaxy (Aharonian et al. 2006a) performed 
with the High Energy Stereoscopic System Cherenkov 
Telescope Array (H.E.S.S.; Hinton 2004). 

Though many galactic VHE sources are recognized to 
be supernova remnants, pulsars and pulsar wind nebulae, 
a small number of them still lack a clear identification. 
HESS J1632-478 is one of them, even if tentatives to ex- 



plain its nature have already been presented (Aharonian 
eral]|2006a| |Walter||2007D . 

A detailed VHE spectral and positional analysis of 



HESS J1632-478 has been reported in Aharonian et al. 
( |2006a[ ), based on 4.5 hours of H.E.S.S. observations. The 
source best fit position is RA=16:32:09 DEC=-47:49:12 
(J2000), placing it in the direction of the near 3 kpc arm 
tangent in the galactic plane. 

The TeV source is extended with a semi-major axis of 
(12 ± 3)^ forming an angle of 21 ± 13° with respect to 



the positive galactic longitude axis. The VHE spectrum, 
between 0.2 and 4.5 TeV, can be fitted with a powerlaw 
model, yielding a photon index T = 2.12 ± 0.20 and a flux 
above 200 GeV of (28.7 ± 5.3) x 10"^^ ^^ ^^-2 g-i 

To understand HESS J1632-478, we collected multi- 
wavelength data from Fermi, XMM-Newton, Spitzer and 
Molonglo, to construct its spectral energy distribution and 
discuss the emission mechanisms. 



2. X-ray Counterpart 

2.1. XMM-Newton observation and data analysis 



The XMM-Newton observatory ( [Jansen et al.||200l| ) in- 
cludes three 1500 cm^ X-ray telescopes each with a 
European Photon Imaging Camera (EPIC) at the focus. 
Two of the EPIC imaging spectrometers use MOS CCDs 
( Turner et al.|20Ql ) and one uses pn CCDs ( jStriider et al. 
2OOTF 



Send offprint requests to: Matteo.Balbo@unige.ch 



In the period from August to September 2008, XMM- 
Newton performed 9 observations of a field centered close 
to HESS J1632-478, collecting data for a total of 92 ks. 
The EPIC observations, listed in table [l] used the thin 
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optical blocking filter. The EPIC pn camera was operated 
in large-window mode. 

Two other XMM-Newton observations of the same 
field were collected in 2003 and 2004. As these observa- 
tions used the medium optical blocking filter and did not 
use the central MOS CCD and the 2004 one was heav- 
ily affected by an outburst of 4U 1630-472, we did not 
combine these observations with those of 2008. 

Table 1. Observation identifiers and date for all 
XMM-Newton observations of 2008 including the source 
HESS J1632-478. The actual and filtered exposure times 
in seconds are given for the two EPIC MOS and for the 
pn camera respectively. 



Obs. ID 



Date 



MOS 

act. s 



MOS 
filt. s 



pn 
act. s 



pn 
filt. s 



556140101 
556140201 
556140301 
556140401 
556140501 
556140601 
556140701 
556140801 
556141001 



14.08.08 
16.08.08 
18.08.08 
20.08.08 
21.08.08 
22.08.08 
24.08.08 
26.08.08 
17.09.08 



10621 


10621 


9834 


9620 


7500 


8034 


8619 


7500 


7034 


12829 


12821 


11234 


4951 


4951 


2239 


14027 


14027 


12398 


10621 


10621 


9834 


11219 


11219 


9634 


9118 


8200 


7534 



5800 
5900 
4500 

11000 
1900 

12398 
6700 
9500 
4500 



Standard data reduction procedures were applied to 
each observation using the XMM Science Analysis 
Software (SAS Version 9.0.0^^ Source positions were 
derived using the SAS task edetect_chain with the EPIC 
MOS and pn data. Only well calibrated single pixel events 
were selected for the pn CCD and single and double 
events for the MOS CCDs. Known hot, or flickering, pix- 
els and electronic noise were rejected. Events were further 
screened with the conservative flag=0 criteria and se- 
lected for the energy ranges (0.2 — 12 keV for MOS and 
0.2 — 15 keV for pn). Periods with enhanced background 

- soft proton flares (EPIC pn count rate above 10 keV 
larger than 0.7 to 4.5 ct/s, depending on the observation) 

- were disregarded in the analysis, resulting in total fil- 
tered exposures of 87 ks and 62 ks for the MOS and pn 
cameras respectively. 

Events were finally selected in the source and back- 
ground regions. Circular regions have been chosen for 
point sources. The regions used for the likely counterpart 
to HESS J1632-478 are described in section l231 

2.2. Serendipitous sources 

In the nine XMM-Newton observations of 2008, three 
point-like sources are always clearly detected: IGR J 16320- 
4751, AX J1632.8-4746 and XMMU J163219.9-474731, by 
order of decreasing X-ray flux. The EPIC derived positions 
of these sources are listed in table [21 



Table 2. Source identifier, name, positions and positional 
statistical error (A) derived from XMM-Newton data. 
The systematic error on the XMM-Newton attitude re- 
construction is of the order of 1.5^'^. 



ID 


Name 


RA 


DEC 


A 


A 


XMMU J163219. 9-474731 


16:32:19.9 


-47:47:31 


0.19^^ 


B 


AX J1632.8-4746 


16:32:48.2 


-47:45:04 


0.13^' 


C 


IGR J16320-4751 


16:32:01.9 


-47:52:29 


0.87^' 


Ee 


HESS J1632-478 


16:32:07.8 


-47:49:09 


1" 


Ep 


HESS J1632-478 


16:32:08.8 


-47:49:01 


0.5'' 



The XMM-Newton data of the high-mass X-ray binary 
IGR J16320-4751 will be presented elsewhere. In these 
data, IGR J16320-4751 features column densities in the 
range (8 — 10) x 10^^ cm~^. 

AX J1632.8-4746 is detected by XMM-Newton as a 
constant source with (14.5 ± 0.3) x 10~^ counts s~^. Its 
spectrum can be well reproduced with the absorbed emis- 



sion of a hot diffuse gas (mekal , Kaastra 1992 ) yielding 
Nh = (6.2 ± 0.2) X lO^^cm-^, kT = 2.5 ± 0.12 keV and 

F2-iokeV = (3.25 ±0.1) X 10-12 erg cm-^ s-^ The ASCA 



data ('Sugizaki et al. 2001) were described with an ab- 



sorbed power-law resulting in a folded model comparable 
to that derived with XMM-Newton. 

The light-curve of XMMU J163219.9-474731 does not 
show any variation. Its spectrum is also well reproduced 
with an absorbed mekal model with Nh = (5.5±0.3) x 10^^ 
cm-2, kT = 2.2 ± 0.1 keV and Fs-iokeV = (l-l ± 0.1) x 
10-12 erg cm-2 s~^. 

The spectra of both AX J1632. 8-4746 and 
XMMU J163219.9-474731 are consistent with these 
of massive stars. If these sources are at a distance > 2 
kpc, their X-ray luminosity (< 10^^ er g/s) is consistent 
with the expectations (10^^-^^ erg/s, Cassinelli et al. 



1981 ). We find possible infrared point-like counterparts to 



these two XMM sources in the Spitzer surveys (Benjamin 



et al. 



2003 Carey et al. 2009) at distances of 3. 3'' and 



2.0^'respectively. 



2.3. HESS J1632-478 Counterpart 

None of the point-like X-ray sources discussed above is a 
likely counterpart of the extended source detected in the 
HESS survey, therefore we merged together all the obser- 
vations to obtain a better sensitivity. To construct mosaic 



images of the MOS and pn data, described in Sect. 2.1 



we used the emosaic command. The mosaic of the count 
images has been divided by the mosaic of the individual 
exposure maps calculated with the eexmap command. 

The resulting MOS and PN mosaic images are pre- 
sented in Fig. [l] An extended source is detected in the 
two mosaic images (source E). The source profile features 
an unresolved point-like component (source Ep in table [2|, 



http://xmm.esac.esa.int/sas/current/ 



^ http://xmm.vilspa.esa.es/docs/documents/CAL-TN- 
0018.pdf 
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with a width (a = 6.8 ± 0.6)'' consistent with the XMM- 
Newton point spread function, in addition to an extended 
component (source Ee in table [2]), close to the detection 
limit. The point and extended sources have significances 
of 15 and 18 a respectively. The extended source could be 
represented by a two dimensional gaussian profile with la 
semi-axes of (32.0 ± 1.2)'' and (15.3 ± 0.7)". The angle of 
the major axis of the extended source is (55 ±2)° counted 
ant i- clockwise from the north direction. 

To extract the X-ray spectra of the point-like and ex- 
tended source, we first defined the source (Fig. [l] bottom) 
and background (Fig.Illtop and middle) extraction regions 
in the mosaic images and selected the events in these re- 
gions for each individual observation. The elliptical region 
used to select the extended source excludes a circle around 
the point-like source (Fig. [l] bottom) . The ancillary spec- 
tral response file for the extended source was generated 
with the extendedsource=YES option of the arf gen tool. 

For the MOS CCDs, the background extraction re- 
gion was chosen in order to minimize the possible con- 
tributions from IGR J16320-4751 and XMMU J163219.9- 
474731 (Fig. IT] top). A different background extraction 
region was chosen for the PN camera (Fig. IT] middle) to 
avoid falling into the gap between the CCD chips when 
extracting the background in individual observations. 

We applied these extraction regions to each observa- 
tion separately. Source and background spectra and re- 
sponse matrices were extracted for each observation. As 
the source is at the same position in the detector plane 
for all the 2008 observations, the resulting spectra were 
combined together using addspec (without error propaga- 
tion, i.e. using Poisson statistics) to obtain a single merged 
spectrum for each camera. 



Table 3. Best fit parameters for an absorbed power-law 
model fitted simultaneously to the three EPIC spectra on 
the point-like and extended sources. Reported errors cor- 
respond to 1(7 uncertainties (68%). The 2-10 keV fiuxes 
are corrected for the absorption. 

Parameter I Point I Extended I Unit 



r 

>-10keV 

..2 




10^ 



10 



-13 —2 - 

erg cm s 



The resulting spectra have been analyzed with XSPEC 
(version ll)|j We fitted simultaneously the merged EPIC 
spectra obtained for the three cameras with an absorbed 
power-law model (Fig. |2|. The resulting xl (table [3| are 
large, probably because the count rate uncertainties re- 
lated to the background subtraction are underestimated. 




Figure 1. The upper and middle panels show XMM- 
Newton EPIC mosaic images in the energy band 0.2-12 
keV for the MOS and 0.2-15 keV for pn CCDs, respec- 
tively. A smoothing algorithm was applied. The height 
of these images is 10 arcmin. Ellipses indicate the source 
(continuous) and background (dashed) extraction regions 
for the extended source. The enlarged bottom image ex- 
tracted from the MOSl camera, shows the extraction 
regions for the point-like (black) and extended (white) 
sources respectively. Events in the black circle are not in- 
cluded in the spectrum of the extended source. The dashed 
white circle shows the region used to select the "back- 
ground" events of the point-like source. 

The absorbed and unabsorbed fiuxes of the extended 



http://heasarc.nasa.gov/docs/xanadu/xspec/xspecll/ 



source are F2- 



lOkeV 



10 



-13 



1.4x10 ^"^ erg cm ^s ^ and 4.3 x 



erg cm ^ s ^, respectively. 
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Figure 2. Top: Best power-law fit of the spectra of the ex- 
tended X-ray source obtained from the three EPIC cam- 
eras. Bottom: 68%, 90% and 99% confidence contours for 
the N// — F parameters. 

We also attempted to fit the point-like source with an 
absorbed black-body model. The derived temperature and 
emitting radius are 1 keV and 22-43 m (for a distance of 
3 kpc). These parameters are not meaningful as a power- 
law model fits the data as well and is probably a better 
representation of the emission at the termination shock. 



3. Radio Counterpart 

The MGPS-2 is a high-resolution and large-scale survey of 
the galactic plane carried with the Molonglo Observatory 
Synthesis Telescope ( MOST) at a frequency of 843 MHz 
(Murphy et al. 2QQ7[). Close to the position of the ex- 



tended X-ray source, the root-mean-square sensitivity is 
1-2 mJy/beam and the beam size has a la semi-axis of 
20^' (obtained by fitting a nearby point-like source). 

Figure [3] displays the MGPS-2 image of the region 
around HESS J1632-478. The extended XMM-Newton 
source lies in a large negative area, caused by artifacts 
of the image reconstruction. An extended diffuse radio 
source, with a size corresponding very well to the XMM- 





1 



Figures. Image (15' x 11') from the second epoch 
Molonglo Galactic Plane Survey (MGPS-2) at 843 MHz. 
Colors are in square root scale between —6 (black) and 18 
(yellow) mJy/beam. The white ellipse shows the region 
used to extract the XMM-Newton detected diffuse source. 



Newton diffuse counterpart (the white ellipse in Fig. [3| is 
clearly observed. By comparing this source with the sur- 
rounding region, a positive excess of the order of ~ 16 
mJy/beam, can be extracted, whereas the distribution of 
the counts in the 4°x 4° image has a la width smaller 
than 3 mJy/beam. This radio excess, with a la semi-axes 
of 35'' X 26", corresponds well to the position and size of 
the extended X-ray source. Integrating this excess yields 
a total fiux density of the order of 19 mJy. A more conser- 
vative estimate of the fiux density, based on the detection 
threshold of the MGPS-2 survey yields to an upper limit 
of 25 mJy. 

We also searched the Parkes 2.4 GHz survey of the 



southern galactic plane (Duncan et al. 1995) performed 



with the Parkes radio telescope for a possible counterpart. 
As the resolution of the Parkes image is only 10.4' and the 
root-mean-square noise is approximately 12 mJy/beam, 
we could only extract an upper limit of ~ 100 mJy at 2.4 
GHz for the extension of the diffuse X-ray source. 

4. Searching the Infrared band using Spitzer 

Moving towards the infrared band we analyzed the 
GLIMPSE and MIPSGAL surveys from the Spitzer Space 
Telescope. 



4.L GLIMPSE 

The Galactic 
Extraordinaire 
was performed with the IRAC instrument (iFazio et al. 



Legacy Infrared Mid-Plane Survey 
(GLIMPSE, [Benj amin et al.| [2003) 



2004) onboard Spitzer, which can simultaneously measure 



in four wavelengths at 3.6, 4.5, 5.8, and 8 /im. The final 
mosaic images have a resolution of ~ 1.2". The 5a point 
source sensitivities are 0.2, 0.2, 0.4 and 0.4 mJy for the 
four different wavelengths, respectively. 

We did not find any sign for diffuse infrared emission 
corresponding to the X-ray and radio counterparts. To 
determine upper limits we analyzed the infrared pixel fiux 
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Figure 4. An extraction from the MIPSGAL survey at 
24 /im (45' X 30'). The small ellipse shows the position 
of the diffuse X-ray source, the dashed ellipse shows the 



HESS source extension (as in Fig. 11 of Aharonian et al. 
2QQ6a[ ) and the two large dashed circles indicates the 95% 
confidence region for the FERMI detections. The crosses 
indicate (A) XMMU J163219.9-474731, (B) AXJ1632.8- 
4746, (C) IGR J16320-4751, (D) the nearby radio pulsar 
PSR J1632-4757 and (F) the unidentified 24/im circular 
structure. 

distribution within and outside of the X-ray diffuse region, 
focusing on the 3.6 and 4.5 /im images. 

We selected all GLIMPSE pixels inside the X-ray 
source ellipse and subtracted the contribution of the 
brightest point-like sources. The final distribution was fit- 
ted with a gaussian and compared with these obtained 
from four other similarly shaped ellipse regions extracted 
in the neighborhood. The average background flux ob- 
tained from the flux distributions varies among the various 
ellipses indicating a gradient of infrared emission in the 
region. The uncertainty on the average background flux 
determination (< 0.04 MJy/sr) integrated over the ellipse 
area represents the minimum flux density necessary to dis- 
tinguish an additional diffuse emission with respect to the 
other regions. For a 3a detection limit, our upper limit 
becomes ~ 25 mJy. 

4.2. MIPSGAL 



MIPSGAL (Carey et al. 2009) is a survey of the inner 



Galactic Plane using the Multiband Infrared Photometer 
for Spitzer (MIPS, ,Rieke et al.|2004D . The survey field was 
imaged in two passbands, 24 and 70 /im with a resolution 
of 6'' and 18'', and an estimated point source sensitivities 
(3cr) of 2 and 75 mJy, respectively. As for the GLIMPSE 
survey, there is no evidence for a diffuse emission corre- 
sponding to the diffuse X-ray source. In addition, at this 
wavelength, the infrared background is strongly inhomo- 
geneous and does not allow to extract a useful upper limit. 
Figure |4] displays the 24 /im MIPSGAL mosaic with 
the various point sources, diffuse X-ray source and error 
circles mentioned in the text. It is worth noting that an 



unidentified circular structure (pointed by the "F" arrow), 
with a radius of ~ 66" ^ visible at 24 /im and not in the 
other bands, lies less than 93'' away from the diffuse X-ray 



source. This object was detected by Flagey et al. (2009) 



but they have not found any counterpart. The nature of 
this object is unclear but its "monochromatic" spectrum 
does not favor non thermal emission and a Supernova rem- 
nant origin. In the 70 and 160 /im there is no evidence for 
any emission coming from the extended X-ray source nor 
from the 24 /im circular structure (source "F"). 

It is also interesting to remark that the extension of 
the TeV source lies in a region of faint emissivity at 24 /im 
(Fig. pi. Perhaps this infrared cavity/edge corresponds to 
dust blown/heated by the progenitor supernova or high- 
energy photons. 

5. The GeV sources in the vicinity 

We looked for possible high energy counterpart in the 



Fermi-LAT first year Catalogue (Abdo et al. 2010), and 



found two unidentified GeV sources in the neighborhood 
of HESS J1632-478 (see Fig. [4|. Both these sources are 
flagged with the letter "c" , indicating that they are to be 
considered as potentially confused with interstellar diffuse 
emission or perhaps spurious detection. Their location, 
flux and spectrum may not be reliable. IFGL J1632.9- 
4802c and IFGL J1632.7-4733c lie 15.5' and respectively 
16.3' far from the extended source detected by XMM- 



Newton. As reported in Abdo et al. (2010), the diffuse 



background model in this region of the galaxy needs to be 
improved and the position of IFGL J1632.7-4733c varies 
with that model. It is therefore very unclear if any of these 
Fermi sources could be a counterpart of HESS J1632-478. 
The HESS error ellipse is parallel to the line joining the 
two Fermi sources, lying close to the extremity of that el- 
lipse. If the two Fermi sources are real and radiate in the 
TeVs, the extension and inclination of HESS J1632-478 
could be significantly affected by confusion. 

6. Discussion 

6.1. Spectral Energy Distribution 

The match in position and size of the radio excess and of 
the extended X-ray source, suggests a non thermal syn- 
chrotron source emitting from the radio to the X-rays. 
The TeV centroid of the H.E.S.S. source lies within the 
extended source detected by XMM-Newton. All other X- 
ray sources are further away and point-like. This positional 
match and the fact that the TeV source is several times 
larger than the X-ray source, suggest that the non thermal 
synchrotron source also emits the VHE emission through 
inverse Compton processes. 

The spectral energy distribution of the source, con- 
structed with the H.E.S.S. and XMM-Newton spectra, in- 
frared upper limits and radio detection or upper limits, is 
featured on Fig. |5J The significant GeV fluxes extracted 
from the Fermi catalogue ( Abdo et al.|2010 ) for the close- 
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by Fermi source FGL J1632.7-4733c are also indicated. 
Together, the XMM- Newton and H.E.S.S. spectra clearly 
indicate the presence of two spectral bumps matching the 
expected synchrotron and inverse Compton emission of a 
Pulsar Wind Nebula. 

If da is the distance of the source in unit of 3 kpc, 
the luminosities of the X-ray extended and point sources 
are ~ 4 x lO^^dg and 2 x lO^^dg erg/s, respectively. 
The bolometric luminosity of the synchrotron and in- 
verse Compton components are lO^^dg and lO^^dg erg/s, 
respectively, assuming the spectrum obtained in Fig. [S] 
Using the empirical relationship of Possenti et al. (2002), 



which presents a high dispersion, the X-ray luminosity is 
as expected for a pulsar with a total spin down power 
of LsD ^ lO^^da"^ erg/s. The ratio of the gamma-ray 
and X-ray fluxes also suggest a pulsar spin- down power of 
3 X 10^^ erg/s and an age of ~ 2 x 10^ years (Mattana 
et al. 2009). With 10% of the current spin down lumi 



nosity radiated at very high energies, HESS J 1632-478 is 
among the oldest and most gamma-ray loud known pul- 
sar wind nebulae CKargaltse v et al.|2007 ), similar as these 
detected in HESS J1825-137 ([Pavlov et al.|2008| ) or HESS 
J1640-465 fLemiere et al.||2009D . 

Although not well defined, the TeV extension of the 
pulsar wind nebula is several times larger than the ex- 
tension of the X-ray nebula, as observed for several aging 
pulsar wind neb ula (HESS J1825-135, HESS J1420-607, 
HESS J1640-465; |Aharonian eTaL]|2006b|c| [Lemiere eFaL 
2009 ) . This is usually explained by the much shorter life- 



time of the electrons emitting synchrotron in the X-ray 
band compared to that of the electrons emitting inverse 
Compton in the TeVs. The very high energy synchrotron 
emitting electrons indeed do not have the time to reach 
the outer parts of the nebula. In the case of HESS J1632- 
478, the TeV extension is larger than the X-ray one, as 
observed in other sources. Note that the value of the ratio 
between the two components extensions can not be well 
constrained due to the low TeV spatial resolution and that 
a better determination could be achieved when more sensi- 
tive/higher resolution Cherenkov instruments will become 
available and possible confusion with the other high en- 
ergy sources detected in the field with Fermi (see Fig. |4| 
could be resolved. 

6.2. Nature of the emission 

An upper limit to the energy density of the synchrotron 
photons in the nebula can be estimated, using the 
peak of the maximum synchrotron emission compatible 
with the observed spectral energy distribution (Fgync < 
lO"-*^^ erg s~^ cm 
of the nebula (0 



^ ) and a lower limit of the angular size 



as 



Us- 



■ sync 



47rcR2 ctan2(0/2) 



< 0.02 eV/cm^ (1) 



where R is the physical size of the nebula. This energy 
density is negligible when compared to that of the CMB 
(UcMB=0.25 eV/cm^) or that of the infrared galactic 
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Figure 5. Spectral energy distribution of HESS J1632- 
478, including the upper limits from PARKES (ma- 
genta) and GLIMPSE (cyan), the probable detection from 
MGPS-2 (green), and the detections by XMM-Newton 
(red) and HESS (purple). The Fermi spectrum of the 
nearby source IFGL 1632.7-4733 is also shown (in blue). 
The continuous line indicates the prediction of the model 
used to represent the emission. At low energy the dotted 
line represents the synchrotron emission from the elec- 
tron distribution described by Eq. (Isl), the dotted-dashed 
and the dashed lines show the CMB and IR dust photons 
components respectively. At high energy, the dotted line 
represents the SSC component, the dotted-dashed and the 
dashed lines show the external Compton emission on the 
CMB and dust photons components respectively. 



background in the galactic plane (Urad ^ few eV/cm^; 
Moskalenko et al.|2006 |) which will dominate the seed pho- 
tons of the Inverse Compton component. 

Assuming that the same electron population is respon- 
sible for the synchrotron and for the inverse Compton 
emission, the ratio of the two components luminosities is 
given by the ratio of the magnetic and radiation energy 
densities 



-^sync 

Lie 



B72/io 
U 



rad 



(2) 



thus the magnetic field in the nebula can be estimated as 



B 



2/ioUrad^^ 



0.5 



Fsync/Fic X 10 /iC. (3) 



As the spectral energy distribution qualitatively indi- 
cates Lgync < Lic/10 (considering the highest energy elec- 
trons, in the synchrotron nebula), the average magnetic 
field in the nebula is B '^ 3 /iC. The spectral energy dis- 
tribution also indicates that the cutoff of the synchrotron 
emission is in the range 10^^~^^ Hz. The energy of the 
electrons emitting at that maximum is in the range 



27r mpC u. 



sync 



eB 



0.5 



(1 - 10) X 10^. (4) 
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The maximum of the inverse Compton emission is 
therefore expected at ~ (1 — 10) x 10^^ Hz, in agreement 
with the TeV spectral energy distribution. 

A possible phenomenological description of the elec- 
tron energy distribution is given by the model recently pre- 



sented by Spitkovsky (2008). This author, starting from 



a numerical two-dimensional particle in cell simulation, 
found a spectrum that, at the downstream of the shock 
front, can be described by a relativistic Maxwellian plus 
a power-law cut-off high energy tail. Such a model has 



been recently adopted by Fang & Zhang (2010) to model 



the broad band SED of a sample of PWNs. We adopt the 
same analytical model that reads: 



Table 4. Emission model representing the spectral energy 
distribution of HESS J1632-478 (assuming a distance of 
3kpc). 





Parameter 


Value 




Electron distribution 


7c 


1.6 X 10^ 






7ci 


1.6 X 10^ 






a 


2.4 






N 


7.0 X 10"^ 


e /cm^ 




R 


1.6xl0^« 


cm 


Synchrotron 


B 


1 


/xG 


External Compton 


ttCMB 


0.25 


eV/cm^ 




ttIR 
'Jrad 


2 


eV/cm^ 



n(7) 



no(7) + i^i(^) 



noil) 

fe) 



exp 



7 < 71 



7> 71, 



(5) 



where 7c is the typical Lorentz factor of the electrons at 
the termination shock, a is the spectral index of the high- 
energy tail, with a cut-off Lorentz factor jci- We fix 71 = 
77c according to "Spitkovsky' (^2008") and Ki is normalized 
in order to match the Maxwellian and the high-energy tail. 

We set the Lorentz factor for the particles at the ter- 
mination shock to 7c = 1.6 x 10^, and use the same value 
of a = 2.4 as obtained in 



6.3. Absorption and distance 

The absorbing column densities measured by XMM- 
Newton in the field are much in excess to the values 
'^H al = (1-5 ~ 1-S) ^ 10^^ cm~^ derived from the low 



resolution radio surveys (jKalberla et al.||2005| jDickey & 
Lockman||1990D . 



The absorbing column density on the pulsar wind 
nebula (Nh ~ 11 x 10^^ cm~^) is very similar to that 
observed for IGR J16320-4751 and slightly larger than 

these observed in the directions of AX J1632. 8-4746 and 

22 



match the X-ray data we use 7ci = 1.6 x 10^. The result- 
ing spectral energy distribution, reproduced using a well 



Spitkovsky] ( |2008D . In order to XMMU J163219.9-474731 (Nh^^x 10^^ cm-^). This in- 
dicates an unusually large column density in this field and 
that about half of the absorption towards IGR J 16320- 
4751 ( [Walter et aT||2QQ 6) may not be related to the High- 
Mass X-Ray Binary system. It also indicates that the dis- 
tance to HESS J1632-478 is of the same order as that of 
IGR J 16320-4751, which was approximately estimated as 
3.5 kpc dRahoui et aT]|2008| ). 



tested code (Tramacere 2007 Tramacere et al. |2009, ) is 
reported in Fig. [5] together with the observed data. We 
reproduce both the synchrotron self Compton (SSC) and 
the external Compton emission, assuming a one-zone (post 
terminal shock) homogeneous emitting region. The details 
of the best-fit model parameters are reported in table [4j 
These parameters are only representative and a simpler 
model could be fitted to the data as well. The quality and 
the uncertainties on the high energy data are such that 
deriving accurate information on the electron distribution 
is not possible. 

Although the Fermi data are partly compatible with 
the emission model, it is difficult to say if the soft GeV 
spectrum of IFGL 1632.7-4733c is related to the pul- 
sar wind nebula or should be considered as upper limits. 
Indeed, the H.E.S.S./XMM source lies well outside of the 
95% confidence region of IFGL 1632.7-4733c (Fig. |4| 

With the shape of the electron distribution inferred 
from the external Inverse Compton component, the 
strength of the synchrotron component is fairly well con- 
strained by the X-ray detection and by the radio detec- 
tion/upper limits. 

The total energy in the electron distribution amounts 
to 10^^ erg. This energy is comparable to the product of 
a spin-down pulsar power at birth of 10^^ erg/s and of a 
characteristic decay time of hundreds of years. 

The nebula magnetic field, size, luminosity and in- 
ferred age are in reasonable agreement with the simula- 



tions by Fang fc Zhang (2010). 



7. Conclusions 

We observed the unidentified TeV source HESS J1632- 
478 with XMM-Newton and looked for counterparts in 
the GeV, infrared and radio bands. An extended faint X- 
ray source is detected close to the centroid of the H.E.S.S. 
error ellipse. A radio excess corresponding to the X-ray 
source is found in the Molonglo sky survey. Upper lim- 
its have been derived from Spitzer and Parkes data. The 
GeV image obtained by Fermi shows two close-by sources 
fiagged as confused in the Fermi catalogue, but none of 
them corresponds to the X-ray source, the situation is 
therefore unclear. 

The ffux density emitted by HESS J1632-478 at very 
high energies is at least 20 times larger than observed 
at the other wavebands probed. The source shape and 
spectral energy distribution suggests a pulsar wind neb- 
ula and can be used to successfully constrain a one zone 
model for the post terminal shock region of the pulsar 
wind nebula. The assumed relativistic electron distribu- 
tion is Maxwellian (7 ~ 10^) with a non thermal tail 
extending to 7 ~ 10^. The synchrotron nebula is faint 
because of the low magnetic field (3 /iG). 
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